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Abstract 

We report on a direct search for a three-body decay of the orthopositronium into 
a photon and two penetrating particles, o-Ps — > 7 + Xi + X2. The existence of this 
decay could explain the discrepancy between the measured and the predicted values 
of the orthopositronium decay rate. From the analysis of the collected data a single 
candidate event is found, consistent with the expected background. This allows to 
set an upper limit on the branching ratio £>(o-Ps — > 7 + Xi + X2)< 4.4 x 10~ 5 (at 
the 90% confidence level), for the photon energy in the range from 40 keV < £" 7 < 
400 keV and for mass values in the kinematical range < mxj + ^x 2 — 900 keV. 
This result unambiguously excludes the o-Ps — > 7 + Xi + X2 decay mode as the 
origin of the discrepancy. 
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1 Introduction 



Positronium (Ps), the positron-electron bound state, is the lightest known 
atomic system. It is bound and annihilates through the electromagnetic in- 
teraction alone. At the current level of experimental and theoretical precision 
this is the only interaction present in this system (see e.g. [1]). This feature 
has made the positronium an ideal system for testing the accuracy of the QED 
calculations for bound states, in particular for the triplet (l 3 ^) state of Ps, 
called orthopositronium (o-Ps). Due to the odd-parity under C-transformation 



the o-Ps decays predominantly into three photons. As compared with the sin- 
glet (l 1 ^) state (parapositronium) decay rate, the o-Ps decay rate, due to the 
phase-space and additional a suppression factors, is about 10 3 times smaller. 
Thus, it is more sensitive to potential admixtures of new interactions, which 
are not accommodated in the Standard Model. 

The study of the o-Ps system has a long history [2]. However, in spite of the 
substantial efforts devoted to the theoretical and experimental determination 
of the o-Ps properties, there is a long-standing puzzle: the o-Ps decay rate in 
vacuum measured by the Ann Arbor group, Y exp = 7.0482 ± 0.0016 [is^ 1 [3], 
has a ~ 5<t discrepancy with respect to the predicted value Y = 7.03830 ± 
0.00007 /is -1 [4] (see also [5]). This discrepancy has been recently confirmed 
by more precise calculations of Adkins et al.[6], including corrections of the 
order a 2 . 

The result of the recent Tokyo measurements of o-Ps decay rate in low density 
SiC>2 powder corrected for matter effects [7] agrees, within the errors, with the 
theoretical value of [6]. However, the method to extract the value of the o-Ps 
decay rate in vacuum from the measurement in matter is still under discussion 
[8]. Thus, it is difficult to disagree with the Adkins et al.[6] statement that : 
...no conclusions can be drawn until the experimental situation is clarified. 

Various exotic o-Ps decay modes have been investigated, with the hope that a 
relative contribution to the o-Ps decay rate at the level of Ar = (Y exp —Y)/Y ~ 
10 -3 would solve the discrepancy^ (for review, see e.g. [8,9]). 

Invisible decays of o-Ps, such as, e.g., o-Ps — > uV, are excluded [10,11]. How- 
ever, there is still a possible explanation of both the discrepancy and the Tokyo 
results. This is based on the existence of the o-Ps — > invisible particles decay 
in vacuum (for more details see Ref. [12]-[14]). 

Visible exotic decays of o-Ps (i.e. decays accompanied by at least one photon 
in the final state) have been experimentally searched for: o-Ps — > 7 + X, 
o-Ps — > 77X and o-Ps — > V7, where X is a new light particle and N = 2,4, ... 
These decay modes have definitely been excluded [15]-[21] as the cause of the 
o-Ps lifetime discrepancy. 

The possibility of a new exotic three-body decay o-Ps — > 7 + Xi + X 2 into a 
photon and a pair of new light weakly interacting particles had not yet been 
considered. The photon energy spectrum from this decay mode has no peak, 
differently from the two-body decay o-Ps — > 7 + X. Thus, the sensitivity of 
previous searches based on the conventional peak-hunting technique is not 
sufficient to exclude this decay mode as the source of the discrepancy. Cos- 

1 Note that the recent result [7] still allows exotic contribution to the o-Ps at a 
(2a) level of Ar ~ 7 x 10" 4 . 
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mological arguments place an indirect stringent constraint on the branching 
ratio of this decay mode [22]. However, one may argue that these arguments 
are depending on many assumptions. 

The purpose of this experiment is to perform a direct search for the decay 
o-Ps — > 7 + Xi + X 2 with the sensitivity needed to exclude this decay unam- 
biguously. 



2 Experimental technique 

The schematic illustration of the detector setup used in the experiment is 
shown in Figure 1. Positrons from the 22 Na source with an activity of 3.6 kBq 
are stopped in the Si0 2 aerogel target (density p ~ 0.1 g/cm 3 , average grain 
size 50 - 100 A) where a fraction of them produces orthopositronium. The 
source is prepared by sealing a drop of the 22 Na solution between two 5 /im 
mylar foils. The source is sandwiched between two 120 /xm thick scintillators 
fabricated by squeezing 1 mm thick scintillator fibers. The light produced by 
positrons crossing one of the scintillator fibers is delivered by the fibers to a 
pair of photomultipliers Philips XP2020 (PMT1 and PMT2 in Figure 1). The 
coincidence of the signals from PMT1 and PMT2 is used to tag the positron 
emission (e + trigger) and to define the time to of positronium formation in 
the target. 

The photons produced by the positronium annihilation are detected by a 4tt 
crystal calorimeter, schematically shown in Figure 1. The detector is com- 
posed of an inner and an outer ring of 8 and 14 BGO crystals, respectively, 
surrounding the target region. Two additional BGO crystals serve as endcaps. 
Each crystal has a hexagonal cross-section with an inner diameter of 55 mm 
and a length of 200 mm. For a more detailed description of the BGO crystals 
see [23]. One of the endcap counters, hereafter called trigger counter, is also 
used for the measurement of the annihilation time of the positronium relative 
to t . When one of the annihilation photons is detected by the trigger counter, 
the recoil photons are detected by the other BGO counters, hereafter called 
the VETO detector. The detector is calibrated and monitored internally using 
the 511 keV annihilation photons and the 1.27 MeV photon emitted by the 
22 Na source in association with the positron emission. Variations of the energy 
scale are within < 1% and are corrected on the basis of an internal calibration 
procedure. 

The particles Xi and X 2 are assumed to be weakly interacting penetrating 
particles. Thus the experimental signature of the o-Ps — > 7 + Xi + X 2 decay is 
the presence of energy deposition in the trigger counter, within a time interval 
consistent with the delayed annihilation of the o-Ps in the target, and no 
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energy deposition in the VETO detector. 

In order to decrease the contribution of two photon events from collisional 
quenching of o-Ps, the Si0 2 target is dehydrated in a vacuum of 10" 2 Torr 
at a temperature about 200°C during two hours before installation inside the 
detector. Furthermore, during the data taking period high purity dry nitrogen 
is flowing through the target. This procedure increases the lifetime of o-Ps in 
the prepared sample of the Si02 aerogel from ~ 70 ns to 132 ns giving a two 
photon suppression factor of ~ 10. 




Fig. 1. Schematic illustration of the experimental setup: a) front view, b) top view. 

The signal of positron emission is also used to open a 3 /is gate for complete 
recording of the signals from the BGO counters and the scintillator PMTs. 
An event is recorded if an energy deposition E TRIG > 40 keV is detected in 
the trigger counter within the 3 /is gate. 

A CAMAC-VME system interfaced to a personal computer is used for data 
acquisition. For each event the following quantities are recorded: 

• the amplitudes, A\ and A 2 , of the pulses from PMTl and PMT2 and the 
time interval Ati2 between them; 

• the energy deposition Etrig in the trigger counter and the time interval 
At e + 7 between the trigger counter pulse and the e + trigger; 

• the pulses from each of the 23 BGO crystals for the measurement of the 
total energy deposition in the VETO detector. 



3 Results 

The search for the o-Ps — > 7 + Xi + X 2 decay described in this paper uses a 
data sample of 2 x 10 7 recorded events. Candidate events for o-Ps decays, with 
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Fig. 2. a) Distribution of the energy deposited by the positrons in the scintillator 
fibers; b) Time difference between the two photomultipliers pulses from the scintil- 
lator fibers; c) Time difference between the scintillator fiber pulses and the trigger 
counter; d) 1.27 MeV photon energy spectrum in the BGO crystals. The dashed 
areas on the plots illustrate the corresponding cuts used for the event selection. See 
Section 3 for details. 

at least one photon in the final state, are identified by the following selection 
criteria, illustrated in Figure 2: 



• the PMT1 and PMT2 pulses from the scintillator fibers in the energy range 
60 keV < Ai t 2 < 400 keV. This cut is applied to reject signals of fake 
positrons from accidental coincidences due to PMT noise; 

• the time difference between PMT1 and PMT2 pulses \At 12 \ < 3.8 ns; 

• the time difference between the PMT pulses from the scintillator fibers and 
the trigger counter signal is required to be in the range 160 < At e + 7 < 
800 ns. The lower cut is used to eliminate background from the tail of 
prompt positron annihilation, shown in Figure 2c). The upper cut is chosen 
to eliminate a region where the background from accidental coincidences 
becomes dominant. 

• the presence in any single BGO crystal, except the two endcaps, of an 
energy deposition in the range 1100 keV < £1.27 < 1500 keV. This criterion 
is applied to select events with a 1.27 MeV photon, emitted by the source 
at the same time as the positron. 
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After imposing the above requirements 338'786 candidates events are found. 
The distribution of these events in the scatter plot Eveto vs -Etr/g is shown in 
Figure 3a). Here Eveto is defined as the sum over all BGO crystals minus the 
energy deposited in the trigger crystal and the energy deposited by the 1.27 
MeV photon: Eveto = — E T rig — -Ei.27- As expected for positronium de- 
cay into photons, the events accumulate around the line EvETO+ETRiG=^'m-e, 
m e =511 keV being the electron or positron mass. 

In Figure 3b) the region of the scatter plot for E yE TO< 35 keV is shown. The 
box for Eveto < 20 keV and E TRIG < 400 keV defines the signal region. The 
cut on the trigger energy Etrig < 400 keV is applied to reject events where 
two photons from the o-Ps — > 37 decay deposit energy in the trigger counter 
and to eliminate background from o-Ps annihilation into two photons due to 
collisional quenching. 

The cut Eveto < 20 keV is chosen according to the measured width of 
the zero-energy peak in the VETO detector. For this measurement, prompt 
positron annihilation events (At e + 7 < 160 ns) are selected with one 511 keV 
photon in the trigger counter. The spectrum measured in the VETO counter 
for zero-energy signal, when the second photon escapes detection, shows that 
95% of the zero-energy signal is collected in the region Eveto < 20 keV. The 
width of the zero-energy signal spectrum is determined mostly by the over- 
lap of close in time events, while the contribution from the ADC pedestals 
fluctuation is found to be negligible. 

The distribution of the VETO energy for Eveto < 500 keV is shown in 
Figure 3c). The distribution of the events with Eveto > 20 keV can be 
extrapolated into the signal region to evalute the background contribution 
in this region. A fit shown in Figure 3c) results in a background estimate of 
Nbckg = 1.6 ± 0.8 events, where the error is evaluated from the uncertainty 
related to the extrapolation procedure itself. 

As shown in Figure 3b) and 3c) one event is found in the signal region. This 
is consistent with the background evaluation described above. Hence, no evi- 
dence for the decay o-Ps — > 7 + Xi + X 2 is found. 

This result allows us to set an upper limit on £?(o-Ps — > 7 + X x + X 2 ) from 
the 90% confidence level (CL) upper limit on the expected number of signal 
events, +Xl+X2 . Because of the uncertainty on the background estimate 

we have chosen conservatively not to subtract the background. Using Poisson 
statistics [24], for 1 event observed and background event expected, the limit 
is C-Ps^ 7 +x 1+ x 2 =3-8 events. 

Given the measured number iV _p s of o-Ps decays in the target, the numbers 
A ir -Ps^7+Xi+x 2 an d No_p s ^s-y of o-Ps — > 7 + Xi + X 2 and o-Ps — > 37 decays , 
respectively, detected in the trigger counter, are 



6 



> 1400 
^ 1200 
o 1000 
m 800 
uf 600 
400 
200 





> 



LU 

> 

LU 



200 



400 



600 



35 
30 
25 
20 
15 
10 
5 






b) 

• • 
* t 


: , , , J i , , 


• 

* . 

• . • .; 



200 



400 



600 



-TRIG 



, keV 



E TRIG> keV 




Fig. 3. Distributions for the selected o-Ps events: a) distribution of the VETO energy 
Eyeto vs the energy Et rig ^ n the trigger counter; b) close-up of the scatter plot for 
Eyeto < 35 keV; the rectangle represents the signal region for o-Ps — > 7 + X1 + X2 
events; c) distribution of Eyeto f° r events in the region 40 keV < -Etr/g < 400 
keV. The line represent the result of a fit of the Eyeto distribution, allowing to 
estimate the background contribution in the signal region. One event is observed in 
the signal region. 

N^ Ps BR{o-Ps -> 7 + X 1 + X 2 )e l7 , (1) 
iVo_P S _»3 7 = ^ _p s e3 7 , (2) 

resulting in 

Br (o-Ps -> 7 + Xx + X 2 ) < ^ o-Ps^ 7+ x 1+ x 2 (3) 

ei 7 A/ _p s ^3 7 . 

Here e i7 and e 37 are the efSciencies of the trigger counter for o-Ps — > 7 + 
Xi + X 2 and for o-Ps — ► 37 decays, respectively. These efficiencies are eval- 
uated with a Monte Carlo simulation of the detector response to both decay 
modes. For the o-Ps — ► 7 + X x + X 2 decay, phase space distribution of the 
photon and of the two particles Xi and X 2 are considered. It is found that 
for the kinematically allowed region of masses < mx 1 + rrix 2 < 900 keV the 
efficiency ratio varies in the range 3.0 < e 37 /ei 7 < 3.7. The lower value corre- 
sponds to ei 7 calculated for m Xl + m X2 =0 keV. The upper value corresponds 
to mxj +mx 2 =900 keV. The efficiency for the o-Ps — > 7 + X1 + X 2 is relatively 
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high over this kinematically allowed region owing to our low trigger threshold 
of 40 keV. For the determination of the limit we take conservatively the value 
£37/ e i7 — 3.7. 

The number of o-Ps decays in the target is measured from the decay curve 
of Figure 2c by fitting the distribution to the function A ■ exp(—t / t _ Ps ) + B 
(B is the accidental background) starting from the time t = 160 ns when 
o-Ps is completely thermalized in the target. Comparing the measured lifetime 
To-Ps = 132.5 ±3.2 ns with the lifetime in vacuum ( 141.9 ns) the probability of 
o-Ps quenching in the target is found to be 6.6%. Correcting for this efficiency 
factor the total number of detected o-Ps is determined to be 3.2 xlO 5 . 

Finally, our 90 % CL limit on I? (o-Ps — > 7 + X x + X 2 ) for the photon energy 
range 40 keV < E 1 < 400 keV and masses m Xl + m X2 < 900 keV is 



5(o-Ps -> 7 + Xi + X 2 ) < 4.4 x 10~ 5 (4) 

This limit is more than 20 times smaller than the value needed to explain 
the discrepancy in o-Ps decay rate. Thus, the o-Ps — > 7 + Xi + X 2 decay is 
definitely excluded as a possible origin of the discrepancy. 
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